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Analytical solutions of and experimental results on the oxi- 
dation kinetics of carbon in a pore are presented. Reaction 
rate, reaction sequence, oxidant partial pressure, total sys- 
tem pressure, pore/crack dimensions, and temperature are 
analyzed with respect to the influence of each on an overall 
linear-parabolic rate relationship. Direct measurement of 
carbon recession is performed using two microcomposite 
model systems oxidized in the temperature range of 700° to 
1200°C, and for times to 35 h. Experimental results are 
evaluated using the derived analytical solutions. Impli- 
cations on the oxidation resistance of continuous-flber- 
reinforced ceramic-matrix composites containing a carbon 
constituent are discussed. 

I. Introduction 

T he kinetics of carbon oxidation are critical for a broad range 
of materials systems, processes, and technologies. And, for 
the specific case of carbon oxidizing and receding in a nonreac- 
tive matrix, the kinetics can be either the Achilles’ heel or an 
enabling step. For example, in a continuous-fiber-reinforced 
ceramic -matrix composite (CFCC), carbon may be present as 
the reinforcing fiber or as a thin film at the fiber-matrix inter- 
face. The interfacial layer provides the critically important 
debond layer which imparts to the composite superior mechani- 
cal properties — most importantly, fracture toughness. The loss 
of the carbon in these materials by oxidation severely reduces 
the service life and is potentially catastrophic. On the other 
hand, the complete removal of the continuous carbon phase 
constituent from a matrix may be desirable in material systems 
where, by design, the loss of “fugitive” carbon forms a benefi- 
cial void. 1 Thus, an understanding of the key parameters con- 
trolling the oxidation kinetics of a continuous carbon phase 
constituent is useful for the design and life prediction of several 
material systems. 

The objective of this work was to answer some of the key 
oxidation questions associated with the technological exploita- 
tion and experimental evaluation of materials with a continuous 
carbon phase. Specifically, can a fiber-reinforced composite be 
designed with a carbon interphase thin enough to slow oxida- 
tion to a technologically acceptable rate? Why are the observed 
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oxidation kinetics in composites coupled to stress? Do carbon 
oxidation kinetics depend on whether CO or C0 2 is formed? 
What is the influence of convection caused by product gases? 
Finally, what experimental variables must be controlled to 
design an experiment that accurately simulates the oxidation 
behavior for a given application? 

Given its historical and current technological importance, the 
oxidation of bulk carbon, in it^ many varied forms, has been 
widely studied. 2 ^ Surface reaction rates are known to be depen- 
dent on a variety of factor:/ such as chemical composition, 
density, porosity and presence of surface defects, crystallinity, 
and crystallographic orientation. In all cases, the rates are 
rapid — in the range of 10“ 4 to 10“ 1 g/(cm 2 -s) which corre- 
sponds to recession rates on the order of mm/h to mm/min) for 
1 atm oxygen at 1000°C. 

Studies on gas phase diffusion-limited oxidation kinetics 
have generally focused on one or a combination of the follow- 
ing three situations: (1) gas diffusion through continuous pores 
or cracks in a “protective” coating on a carbon substrate (e.g., 
coated carbon/carbon composites); 5 8 (2) pore/crack sealing by 
a condensed matrix-oxidation reaction product; 7-11 or (3) oxida- 
tion of a thin carbon film in a nonreactive matrix for a specific 
condition of interest. 1213 Of particular relevance to the work that 
follows, prior studies have shown that 800°C is an important 
temperature when considering oxidation of carbon fiber-rein- 
forced composites. At temperatures below 800°C, the rate of 
carbon recession is dominated by the kinetics of the chemical 
reaction at the receding carbon surface. In contrast, at tempera- 
tures greater than 800°C (but low enough that pore/crack seal- 
ing is not an issue), gas phase diffusion assumes an increasingly 
important role in determining kinetics. This paper generalizes 
oxidation of thin carbon films and fibers in nonreactive matrices 
with a formal analysis, explicitly including both chemical reac- 
tion and gas phase diffusion effects. The resulting general 
model for oxidation kinetics provides the basis for evaluating 
other limiting cases. 

In the following section, a comprehensive analytical model 
is developed using reaction rate formalism and kinetic theory 
of gases. The results of two sets of experiments are then inter- 
preted using aspects of this model in Section III. In the experi- 
ments, recession distance is measured directly, using either 
metallography or precise measurements of dimensional 
changes. This is in contrast to most oxidation studies, where 
oxidation rate is measured indirectly, i.e., through degradation 
in mechanical properties, weight changes, or changes in elec- 
trical resistivity. 3,13,14 The advantage of directly measuring car- 
bon recession is that it provides the opportunity to validate the 
assumptions used in derivation of the model. 
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II. Modeling of the Oxidation Process 

The variables affecting the rate of carbon oxidation include 
the kinetics of the oxidation reaction occurring at the carbon 
surface, the characteristic dimension of the carbon (the fiber 
diameter or the coating thickness), the nature of the reaction 
product (either CO or C0 2 ), the molecular weight of the gases, 
the temperature, and the total ambient pressure. One advantage 
of a rigorous analysis is that the sensitivity of the oxidation 
process to each of these variables can be assessed. The analysis 
focuses on obtaining an expression for the carbon recession 
rate. In all cases, it is assumed that the dense and homogeneous 
carbon initially occupies a well-defined volume which is finite 
in the two directions which are perpendicular to the oxidation 
direction and infinite in the direction parallel to the oxidation 
direction, as shown schematically in Fig. 1 . The results of the 
model are equally applicable to the oxidation of a carbon fiber 
(in a carbon-fiber-reinforced CFCC) and a carbon coating (for 
example, the fiber/matrix interface coating in a silicon carbide- 
fiber-reinforced composite). The discussion below is limited to 
the oxidation of a thin annulus of carbon between a nonreactive 
fiber and matrix. For silicon-based CFCCs, oxidation of both 
the fiber and matrix can be neglected at “intermediate” tempera- 
tures (less than ^ 1000°C). This assumption has been used 
because it has been shown that oxidative damage is most severe 
in this intermediate temperature range. 1415 It is also assumed 
that the receding carbon maintains a planar interface. 

The assumption is made that the oxidation of carbon is gov- 
erned by sequential linear-parabolic kinetics, i.e., the recession 
distance, x, is related to the oxidation time, t, by 16 

x 2 /k p + x/k x = t (1) 

where k p is termed the parabolic rate constant (m 2 /s) and k x the 
linear rate constant (m/s). When the recession distance is small, 
the left-hand side of Eq. (1) is dominated by the linear term, 
but as the recession distance increases there will always be a 
transition to parabolic kinetics. It is common to refer to the 
former as reaction control and to the latter as diffusion con- 
trol, although some diffusion-controlled processes give linear 
kinetics and some reaction-controlled processes are nonlinear. 
In the analysis that follows, the two limiting cases are treated 
independently in order to isolate the importance of system 
variables on each. 

(1) Reaction-Controlled Kinetics 

In this case, the reaction at the receding carbon interface is 
assumed to be first order, and either 

0 2 + C <-> C0 2 (2a) 

or 

0 2 + 2C <-» 2CO (2b) 

The rate of the reaction, r (mol/m 2 - s), is given by 

r = Kfii (3) 

where K, is the reaction rate constant (m/s) for the oxidation 



x=4 


Fig. 1 . Schematic of carbon recession during oxidation in a nonreac- 
tive matrix. 


reaction and c t is the concentration of oxygen (mol/m 3 ) at 
x ~ ^ (see Fig. 1). If it is appropriate to consider the diffusive 
supply of oxygen to be infinitely fast, then c € is effectively 
equal to c D , the concentration of oxygen in the atmosphere. c Q 
can be alternately written as xc T > where x is the oxidant partial 
pressure and c T is the total concentration of gas molecules 
the atmosphere (mol/m 3 ). In this case, the carbon reces 
rate is constant (i.e., independent of time) and so-called li 
kinetics are observed, i.e., 

dtydt = k x 
Equivalently, 

£ = (5) 

if no recession has occurred prior to t = 0. The linear rate 
constant can be written as 

K = K r xc T /N (6) 

where N is the molar density of carbon (mol/m 3 ). The reaction 
rate constant is conventionally written using an Arrhenius form, 
and c T can be written using the ideal gas law. Thus, an alternate 
expression for k x , in terms of fundamental quantities, is 

K = (\IN)( X PiRT)k 0 exp (-Q/RT) (7) 

where k 0 is a constant (m/s), Q is the activation energy (J/mol), 
P is the pressure (Pa), R is the gas constant (J/moTK), and T is 
absolute temperature (K). Values for N, k 0 and Q are dependent 
on the nature of the carbon and must be determined experimen- 
tally. There is a wide variation in the reported values for oxida- 
tion rates of bulk carbon. 2,3 Assuming these experimentally 
determined values are constant over the temperature range of 
interest, the functional dependence of temperature and pressure 
on the linear oxidation rate can be predicted from Eq. (7). For 
pressure, the oxidation rate is expected to vary linearly. The 
temperature dependence is readily observed by considering the 
product of the linear rate constant times temperature, or 

k x T exp(~Q/RT) (8) 

When the argument of the exponential function is small 
(e.g., activation energy is much less than the product RT) the 
k x T product approaches a constant value, i.e., inverse linear 
dependence. On the other hand, when the exponential argument 
is large, the reaction rate increases exponentially with increas- 
ing temperature. 

(2) Diffusion-Controlled Kinetics 

In this case, the reaction at the carbon interface is assumed to 
be infinitely fast, but several distinct reaction sequences are still 
possible and affect mass transport. For example, the case of 
simultaneous oxidation of carbon and silicon carbide in CFCCs 
was considered in a recent set of papers. Two cases were ana- 
lyzed: 0 2 diffusion in and C0 2 diffusion out, and 0 2 diffusion 
in and CO diffusion out. 9 It is also possible for CO formed at 
the carbon interface to undergo a secondary oxidation to form 
C0 2 . This case has also been analyzed. 12 The gas concentration 
profiles in a crack/pore for each of the three reaction sequence 
cases are shown in Figs. 2(a), (b), and (c). All three will be 
considered here. It should be noted that although only one 
reaction sequence can be the correct description for a particular 
situation, it is plausible that different sequences will occur 
under different conditions. 

It greatly simplifies the analysis to treat the diffusion coeffi- 
cients of all of the species in the gas phase as identical in 
magnitude. This is not strictly true, but the differences are 
small, as discussed below. 

Gas phase diffusion coefficients can be predicted using 
Chapman-Enskog theory. The expression for the diffusion 
coefficient of gas A, in a stationary atmosphere of B, is 17 

D ab - 5.9543 X 1(T 24 [(1/M a ) +(1/M b )] 1/2 T 3/2 / 

(P& ab 2 ^ab) 


( 9 ) 
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Chemical Reaction 


(a) 0 2 + c ~~ C0 2 


(b) O 2 + 2C 2CO 


Gas Concentration Profile 



(c) 


At 4 

C0 2 + c — -2CO 

At £/f 

0 2 + 2CO 2C0 2 



Fig. 2. Plausible oxidation reaction sequences: (a) 0 2 diffusion in 
and CO z diffusion out, (b) O z diffusion in and CO diffusion out, (c) CO 
formation at the carbon surface with subsequent oxidation to C0 2 . 


fixed by the composition of the atmosphere to be \ c t and the 
concentration at the interface, x = is fixed by the appropriate 
phase equilibria (Eq. (2a)) to be c ? . It is further assumed that 
the quasi steady state exists, i.e., the flux density of oxygen,y ox , 
is not a function of position at any time t. The geometrical 
arrangement is schematically illustrated in Fig. 2(a). The flux 
density of oxygen over the domain 0 < x < £ may be, in 
general, written 18 

7 ox '= ~D dcjdx + v*c m (13) 

where v* is the molar average velocity, which may be defined 
using 

v* = (Xc,v,)/(Xc ( ) (14) 

where v, is the velocity of the /th species with respect to a 
stationary coordinate system. By definition 

ji = c,v, (15) 

so that Eq. (14) can be rewritten, in the general form, as 

v* = X(i,)/c T (16) 


where D is the diffusivity (m 2 /s), T is the absolute temperature 
(K), M is the molecular weight (kg/mol), P is the total pressure 
(Pa), cr is the collision diameter (m), and fl is a tabulated 
integral which is a function of the energy of molecular interac- 
tion parameter, e (typically tabulated as a ratio to Boltzmann’s 
constant, k). For mixtures, the effective collision diameter and 
energy of interaction are given by 

o- AB = (o'a + o- b )/2 (10) 

and 

E A Jk^[(E A !k)(eJk)V /2 (11) 

respectively. Inspection of Eqs. (9), (10), and (11) indicates 
that D AB = D ba . For the cases under consideration here, it is 
convenient to examine the ratio D AB /D AC where A, B, and C 
refer to different gas species. With Eq. (9) it is possible to 
express this ratio as 

D a JD ac = {[M c (M a + M b )]/[M b (M a + M b )]}' /2 

[(ct a + a c )/(a A + cr B )] 2 (fl AC /0 AB ) (12) 

The relevant parameters for the systems of interest are given 
in Table I. Ignoring the slight temperature dependence of 
Cl, and substituting these values into Eq. (12) yields D coo J 
^co,co 2 — 1.28, Do 2 ,col^o-,,co 2 — 1-30, and D COiO JD C0 ^ C0 — 
0.99. The temperature dependence of fl changes these ratios by 
less than 1% over the temperature range of interest (800°- 
1200°C). In other words, the diffusion coefficients vary by less 
than 30%, even in the most extreme cases. 

Conservation of the oxidant proves to be the most convenient 
frame of reference for most of the discussion. It is shown that 
the solution for the most complex case, CO formation with 
subsequent oxidation to C0 2 , is just a combination of the two 
simpler solutions for CO and C0 2 as single reaction products. 
Further generalization to include Knudsen effects (i.e., gas dif- 
fusion limited by dimensions of the channel) is then offered. 

(A) 0 2 Diffusion In and C0 2 Diffusion Out: It is assumed 

that the concentration of the oxidant at the surface, x — 0, is 


where c x is the total concentration of gas molecules, which 
is a constant. For this particular case, where c T is simply the 
sum of the concentrations of oxygen and carbon monoxide. 


Eq. (16) becomes 

v* = O'o* + 7 co 2 )/ 6 't ( 17 > 

Substituting Eq. (17) into Eq. (13) yields 

j 0 , = -D dcjdx + (cJc T )(j ox + j c J (18) 

Conservation of mass requires 

7co 2 = -/ox (19) 

and therefore v* = 0 and Eq. (13) simplifies to Fick’s first law 

jo x = ~D dcjdx (20) 

Assuming quasi steady state yields a linear concentration gradi- 
ent, which, upon integration gives 

7ox - -D(cJ0 - cjm/i (21) 

The recession rate is given by mass balance; i.e., the thickness 


removed is equal to the negative of the flux density of carbon 
away from the interface divided by the molar density of carbon. 


N, (mol/m 3 ) in the condensed state, or, 

d£/dt = JJN (22) 

Substituting Eqs. (19) and (21) into (22) gives a standard form 
for parabolic kinetics 

dUdt = [D (c ox (0) - c ox (0)]/(/V0 (23) 

and integrating gives 

e = [2D(c ox (0)-c ox (Q)/N]r (24) 

which can be written in familiar notation as 

€ = K /2 t' /2 (25) 


Table I. Parameters Used in Calculation of Gas Phase Diffusion Coefficients Using 


Chapman-Enskog Theory* 



M 

a 

e/k 


Vl(e, /*)(£,/*)] 

’Qab 

Gas 

(kg/mol) 

(m) 

(K) 

Mixture 

(K) 

(at 1073 K) 

0 2 

32 X 10- 3 

3.467 X KT 10 

106.7 

o 2 /co 

99 

0.737 

CO 

28 X 10- 3 

3.690 X 10 10 

91.7 

co/co 2 

144 

0.784 

C0 2 

44 X 10“ 3 

3.941 X 10-'° 

195.2 

co 2 /o 2 

134 

0.774 


*Vaiues for 11 are interpolated from tabulated values.’ 7 
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with k p , termed the parabolic rate constant, given by 

K = [2£>(c ox (0) - c m (0)/N] (26) 

(B) 0 2 Diffusion In and CO Diffusion Out: In this case it 
is assumed that reaction (2b) occurs at x — £ and that the CO 
diffuses, without reaction, to the external surface of the pore. 
This is plausible, since it is well known that CO may be main- 
tained as a metastable species for long times. The concentration 
gradients for this case are illustrated in Fig. 2(b). The chemical 
reaction in this case produces two moles of gaseous product, 
CO, for.every mole of gaseous reactant, 0 2 , consumed, i.e., 

Ac = -2/ ox (27) 

which makes the convection term in Eq. (13) nonzero. Substi- 
tuting Eq. (27) into Eq. (16) and in turn into Eq. (13) yields, 
upon rearrangement, 

jox = -KDc t )/(c t + cj\ dcjdx (28) 

Integration of this expression does not yield a linear concentra- 
tion gradient. The assumption of quasi steady state, i.e., j ox is 
not a function of position at any time t, allows the right-hand 


side of Eq. (28) to be set equal to a constant 

A = —[Dc t I(c t + cj] dcjdx (29) 

Since d(c r + cj/dx = dc ox /dx, this may be rearranged to give 
<5 [In (c T + cj]/dx = A' (30) 

and upon integration 

In (c T 4- c ox ) = A'x + B (31) 

or 

c ox = B’ exp (A'x) - Cy (32) 

Given the boundary conditions 

cjx = 0,0 = c ox (0) (33) 

Co*0 = €.0 = cJO (34) 

the constants can be evaluated to give 
Cox = (C T + c ox (0)) 

exp {In [(c T + c m (ZMc r +cj O))](x/0} 

- c T (35) 

Upon substitution into Eq. (28), this yields 

j ox = -(Dc r /Q In [(c T + c ox (£))/(c x + c ox (0))] (36) 

The rate of carbon recession in this case is 

dm = ~j co /N = 2 jJN (37) 

Substituting Eq. (36) and integrating yields 

i 2 = (4 Dc t /N) In [(c T + c ox (0))/(c T + c ox (0)]t (38) 

The kinetics remain parabolic, but the parabolic rate constant, 
k p , is now defined as 

K = (4 Dc r IN) In [c r + c ox (0))/(c T + c ox (m (39) 
or equivalently 


K = (4 Dcy/M) In [(1 + (c ox (0)/c T ))/(l + c ox (0/c T ))] 

(40) 

(C) Formation of CO at Receding Carbon Surface with 
Subsequent Oxidation to C0 2 : In this case it is assumed that 
the atmosphere is pure 0 2 and that CO is formed at the receding 
surface. Following the analysis of Bernstein and Koger, 12 the 
CO is oxidized to form C0 2 at a position £//, where / is a 
constant. The relevant chemical reactions are 

(41) 


and 

2CO + 0 2 -> 2C0 2 at x = £/f (42) 

The concentration gradients of the gaseous species are sche- 
matically illustrated in Fig. 2(c). The previous analysis was 
from the (valid) point of view of continuity of carbon flux. 12 
The analysis given below focuses on the flux of oxidant, in this 
case C0 2 , and is intended to illustrate that this problem is 
a straightforward extension of the two simpler cases given 
in Sections II(2)(A) and II(2)(B). It is convenient to divide 
the problem into two regions: Region I includes the domain 
0 ^ x ^ £//; and Region II is defined by £//< x < Note that 
iff = 1 then Region I covers the entire range 0<xS(, and 
the system is equivalent to that presented in II(2)(A), and if 
f =(X> then the entire spatial domain falls into Region II and the 
case discussed in II(2)(B) applies. 

In Region I, carbon dioxide is the oxidant and its flux is 
described in a manner analogous to Eq. (13), giving 

j C o 2 = ~D dc c Jdx + v*c c02 (43) 

but this again simplifies to Fick’s first law. The reasoning for 
setting v* equal to zero in this case is that although reaction 
(41) consumes 1 mol of gaseous reactant to make 2 mol of 
gaseous product, this “extra” mole is consumed in reaction (42), 
which yields only 2 mol of gaseous products from 3 mol of 
reactants. Effectively, there is a negative convective term at 
x = £ which is balanced by a positive convective term of equal 
magnitude atx = £//. Therefore, quasi steady state again yields 
a linear concentration gradient and 

l 02 = ~D[c c jm - c C02 m/(m (44) 

and since 

dm = JJN (45) 

substitution and integration yields 

e = (iDfmcjm - cjm ( 46 ) 

yielding a parabolic rate constant, k" p , for this case 

K = C 2DflN)[c C02 (i/f ) - c CO2 (0)] (47) 

which involves the unknown 

Turning attention to Region II,/ is given by Eq. (43), but in 
this case v* is nonzero and analogous to Eq. (28). Substitution 
for v* using fluxes, in this case / and / co , and rearranging, 
gives 

jco 2 = ~[(DCy)/(Cy + c c J} dc C0 Jdx (48) 

Applying the assumption of quasi steady state and the following 
boundary conditions 

Coo.U = £//,() = (49) 

<•„,.(* = CD = c„,(i) (50) 

gives 

<4o 2 = (c C o 2 (0 + c t) 

exp {[//(/- 1)] In [(c T + cji))/ 

(Cy + c C02 (W))]Wt- 1)} -c T (51) 

Substitution of Eq. (48) into Eq. (45) and rearrangement yields 

/ C02 = -(DcylO [//(/- 1)] 

In l(Cy + c CO2 (0)/ (C T + cjum (52) 

Now, for the case of this reaction sequence, the recession rate is 
given by 

dm = ~l/2jJN = j c JN 


CO, + C — > 2CO at x = ^ 


(53) 
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The origin of the one-half in Eq. (53) comes from the fact that 
only one of the CO molecules on the right-hand side of reaction 
(41) removes a carbon from the condensed carbon. Substitution 
of (52) in (53) and integration yields 

e = (2 Dc t /N) [//(/ — 1)] 

In [(C T + c C02 m))/(Cr + c C02 m t (54) 

and an alternate form for the parabolic rate constant, k' p 

K = (2 Dc T IN)[fKf- 1)] 

In [(c T + c C02 W))/(c t + c COi m , (55) 

Equating Eqs. (47) and (55) and solving for “/” yields 

/= [c T /(c C02 m) - c C02 (m 

In [(c T + c C02 (ilf ))/(c T + c C02 m + 1 (56) 

and substituting back into either (47) or (55) allows “/” to be 
eliminated from the expression for the parabolic rate constant, 

K = VD/N){c t In [(c T + cjf/f))/ (c T + c C02 m 

+ k CO2 (^//)-c CO2 (0)]} (57) 

(D) Comparison of the Parabolic Rate Constants: At first 

inspection, Eqs. (26), (40), and (57) appear quite different. 
However, referring to the schematics (Figs. 2(a-c)) illustrates 
how these may be simplified and compared. For the case of 0 2 
diffusion in and C0 2 diffusion out, c ox (0) x c t and c ox (£j) ^ 0 
so that 

K ~ 2 X (Dc T /N) (58) 

Similarly, for the case of 0 2 diffusion in and CO diffusion out, 
Cqx(O) *= xc T and cjf) *=» 0, so that 

£ p ^ 4 In [1 + x](Dc t /N) (59) 

Finally, for the case of CO formation at receding carbon 
surface with subsequent oxidation to C0 2 at £//, setting 
c C02 m ~ X^T, c CO2 (0) « 0 and c COj (0 0 yields 

/= 1 + (1/x) In (1 + x) (60) 

and 

*p~ 2[ X + ln(l + X )KDc t /N) (61) 

Equations (58), (59), and (61) differ only in the magnitude of 
the numerical constants which precede the common Dc r /N 
factor. It is interesting to compare their magnitudes and reflect 
on their physical interpretation. 

0 2 diffusion in and C0 2 diffusion out yields the slowest 
kinetics for all values of x- 0 2 diffusion in and CO diffusion out 
predicts kinetics which are more rapid, but only modestly. For 
example, for 1 atm pure 0 2 (x - 1), k' p /k p 1.39. Mass balance 
would suggest that the formation of CO would remove carbon 
twice as fast as the formation of C0 2 , since 2 mol of carbon 
would be removed for each mole of oxygen which diffuses in. 
The net increase, however, is only 39%, because the convection 
caused by the outward flow of CO opposes diffusion of the 
oxidant. This modest increase in oxidation kinetics would be 
further reduced when the additional effect of variation in the 
gas phase diffusion coefficients, Eq. (12), is taken into account. 
For example, if one assumed no gradients (maximum effect), 
the ratio k p /k p becomes nearly unity. 

The formation of CO with subsequent oxidation to C0 2 
yields a prediction of the most rapid kinetics, for 1 atm pure 0 2 . 
Still, the differences are small; C[Jk p *=* 1.69 and k' p /k' p *=» 1.22 
(ignoring diffusion coefficient variations). It may seem counter- 
intuitive for this case to yield the most rapid kinetics, since a 
global mass balance states that only 1 mol of carbon will be 
removed for each mole of oxygen. The physical reason for the 
fast kinetics is that the existence of two regions yields steeper 
concentration gradients, for all species, at a given recession 


distance. This is illustrated by comparing Figs. 2(a), (b), and 
(c). Faster recession occurs because the kinetics of carbon 
recession are directly proportional to the flux of oxidant, which, 
in turn, is directly proportional to the oxygen concentration 
gradient. 

(3) Knudsen Effects 

The calculations carried out in Section 11(2) involved the 
implied assumption that gas phase diffusion was determined by 
intermolecular collisions (Eq. (9)), which is sometimes termed 
“molecular diffusion.” It is plausible that if carbon is present as 
either very small-diameter fibers or very thin interlayers, gas 
phase diffusion will be affected by collisions between mole- 
cules and the walls of the pore left by the receding carbon. The 
limiting case where molecule-wall collisions are much more 
frequent than intermolecular collisions is termed “Knudsen 
diffusion,” and is characterized by a significantly smaller diffu- 
sion coefficient. 

It is perhaps most useful to analyze the transition case, which 
is general, and take limiting cases to recover the expressions 
appropriate for the molecular and Knudsen regimes. An 
existing analysis may be readily extended to yield the parabolic 
rate constant for this case. 17 For the direct formation of C0 2 , 
Section II(2)(A), the Knudsen case yields a parabolic rate con- 
stant of form identical to Eq. (26) except that an effective 
diffusion coefficient, £> eff , must be used. D eff is given by 

1/Deff = 1 ID + 1 ID k (62) 

It is trivial to take the limiting cases. 

Both of the other plausible reaction sequences involve the 
formation of CO at the receding interface. Since the premise of 
the analysis is that intermolecular collisions are rare events, it is 
assumed that the case analyzed in II(2)(B), i.e., CO diffusion 
out without subsequent reaction, applies, and this is the case 
analyzed here. It is also possible to extend the analysis to 
include the case of CO formation with subsequent reaction to 
C0 2 , but we have elected not to, since it is unlikely that gas 
phase equilibrium will be attained over distances of interest 
in CFCC’s. 

An expression has been derived 17 which describes the flux 
density of a species in a binary gas mixture in the transition 
regime, which may be written in the notation used here as 

Ax = -c T l[((l - a cJ/D) + (1/D k )] dcjdx (63) 

where D k is the Knudsen diffusion coefficient which, for a 
cylindrical pore, can be calculated using 17 

D k = 1.534 dj(T/M) (64) 

where d is the diameter of the pore (m), T is absolute tempera- 
ture (K), and M is molecular weight (kg/mol). D k for an annular 
pore of thickness d' (such as would be left by a receding 
interlayer) would be modestly larger, because the average dis- 
tance between collisions would be larger than d. 

The parameter a in Eq. (63) is dependent on the relationship 
between the fluxes of the two species. For the oxidation case 
being analyzed here, 

a = 1 + jJjo* (65) 

and, by Eq. (27), a ~ - 1, so that Eq. (63) is simply 

j ox = -c r /m + cJ/D) + (1 /D k )\ dcjdx (66) 

Substitution in Eq. (37), integrating and rearranging gives 
e = (4 Dc r IN) 

In {[(1 + cJ0)/c t )(DJD) + 1]/ 

[(1 + c 0 M)/c r )(DJD) + l])r (67) 

The generalized expression for k' p is therefore 
k' p = (4 Dc r /N) 

In {[(1 + cJ0)/c T )(DJD ) + 1]/ 

[(1 + Co M/Cr)(DJD ) + 1 ]} 


( 68 ) 
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or equivalently 
k' p = (4DcJN) 

In {[(1 + c ox (0)/c x ) + (D/D k )]/ 

[(1 + c ox (0/c T ) + (D/D k )]} (69) 

The limiting cases are a function of the ratio of the molecular 
to Knudsen diffusion coefficients. The molecular regime is 
obtained when D « D k or D/D k — > 0, and it may be seen by 
inspection that Eq. (40) is recovered from Eq. (69) when this 
substitution is made. 

The other limiting case is most clear when the approxima- 
tions c ox (0) «=» x c t and c ox(£) w 0 are made. With these, 
Eq. (68) may be rewritten as 

K = ( 4 Dc t/N) In {[(1 + X )(DJD) + 1]/ 

KDJD)+1]} (70) 

In the Knudsen regime, DJD is small and k p may be validly 
approximated, using In (1 + x) x, for small x, as 

K ~ (4DcJN){(\ + X )(DJD) - (DJD)} 

~ (4D X cJN)(DJD) = 4 X (D k c T /N) (71) 

Thus, the value of k p in the Knudsen limiting case is differenti- 
ated from that in the molecular case, Eq. (59), by both the 
magnitude of the diffusion coefficient and numerical factor. The 
differences are small. The value of the numerical coefficient in 
the Knudsen case is higher than the molecular case because of 
the absence of intermolecular collisions that preclude convec- 
tion effects. In a pure oxygen atmosphere, i.e., X = 1, the 
numerical coefficient in the Knudsen case is a factor of 1.69 
larger than that in the molecular case, which tends to moderate 
the effect of the difference between D and D k . As x is reduced, 
the values of the numerical coefficients converge. 

Since D k is dependent on the dimensions of the pore, k' p may 
be a function of, for example, interface thickness. It is possible 
to estimate the value of as a function of the characteristic 
dimension of the pore using Eq. (70). For 0 2 , D k (m 2 /s) is 
obtained from Eq. (64). 

D k = S.515T i/2 d (72) 

c T , in mol/m 3 , is calculated from the ideal gas law, 

c T = P/(RT) = 1.203 X 10 ~\P/T) (73) 

N was calculated assuming the density of carbon, which varies 
depending on the degree of crystallinity, to be 1.8 g/cm 3 , or 

N - p/M = 1.5 X 10 5 mol/m 3 (74) 

D, in m 2 /s, can be estimated from Chapman-Enskog theory. A 
binary mixture of CO-0 2 was assumed in the calculation. From 
Eqs. (9), (10), and (11), and the values in Table I, 

D = 1.593 X 10" 4 (r 3/2 /F) (75) 

Combining these expressions with Eq. (70) gives 

k’ p = 6.263 X l0-'°r 1/2 

In {((1 + x)4-392 X l0\Pd/T) + 1)/ 

(4.392 X 10 \Pd/T) + 1)} (76) 

A family of three curves generated using this function, with 

X = 1 (1 atm pure 0 2 ), are shown in Fig. 3. Two variables are 

explored: temperature and pressure. It is evident that tempera- 
ture only weakly affects the magnitude and shape of the curve. 
In the molecular (large d) regime, D is proportional to T m and 
c T is inversely proportional to temperature, so k' p , is directly 
proportional to the square root of temperature. At 1 atm pres- 
sure, the limiting values in the molecular regime are calculated 
to be 1.42 X 10 8 and 1.55 X 10" 8 m 2 /s, for 800° and 1000°C, 
respectively; i.e., the difference is less than 10%. These values 
for k p predict recession of 1 cm in roughly 2 h — 7030 and 6460 



Fig. 3. Effect of variations in temperature and pressure on the para- 
bolic rate constant, k p (assuming 0 2 + 2C — > 2CO). 


s, respectively. The temperature effect in the Knudsen regime is 
also modest, but, interestingly, it is of opposite sign because D k 
has a weaker temperature dependence, T 1/2 , than D. In this 
regime, k p is inversely proportional to the square root of temper- 
ature; at higher temperatures, k p is smaller. 

Secondly, a reduction in the pore opening diameter, d, to 10 
|xm is predicted to produce almost no effect for any of the three 
cases considered. More importantly, the magnitude of Knudsen 
effects is always small; at 1000°C and 1 atm pure 0 2 , when “af” 
is 1 |xm, k p is 80% k p (oo) and when “d” is 0.1 |xm (which is the 
thickness often quoted as the minimum permissible for desir- 
able mechanical properties in CFCCs 19 ) k p (().l)/k' p ( r ~r>) is 1/3. 
Finally, increasing pressure does not change the value of k' p in 
the molecular regime, because the direct proportionality of c T is 
nullified by the inverse proportionality of D. Flowever, at higher 
pressures the already modest Knudsen effects are further sup- 
pressed (at 15 atm and d = 0.1 |xm, k' p and k p (°°) differ by less 
than 2%). 

III. Experimental Procedure 

Oxidation experiments using two model composite systems 
were conducted. Each consists of a carbon cylinder in a nonre- 
active matrix. Schematic cross sections of each are shown in 
Fig. 4. The two are principally distinguished by the type of 
carbon and 2 orders of magnitude difference in diameter. These 
diameter values are indicated on the scale in Fig. 3. It can be 
seen that k p for these two systems should be indistinguishably 
close. Further, either system gives a k p that closely approxi- 
mates the value expected for a 2 p,m thick interphase, (i.e., 
<20% different), provided the ambient pressure is held con- 
stant at 1 atm. Thus, the analytical model allows results from 
model systems to be applied to carbon recession in composites, 
where carbon is used either as a fiber or interphase. 



3.175 mm 33 ^ 

(a) (b) 


Fig. 4. Schematics of model microcomposites used in oxidation 
experiments: (a) carbon rod in A1 2 0 3 , (b) pitch-based carbon in silicon 
carbide sheath (SCS-6 fiber). 
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The first model system was constructed by inserting 3.175 
mm ( 1/8 in.) diameter carbon rods into hollow alumina cylin- 
ders. The rods are commercially available (E. F. Fullum, Sche- 
nectady, NY) and commonly used as electrodes for preparing 
physically vapor deposited carbon films. Both the rods and the 
alumina had initial lengths of nominally 65 mm. Individual 
assemblies were exposed at temperatures in the range of 700° 
to 1200°C and for times to 24 h in a muffle furnace. Following 
exposure, the carbon rod was removed from the alumina sheath 
and the dimensions were measured. The recession distance was 
calculated by halving the difference between the length of the 
alumina tube and that of the carbon rod after exposure. In all 
cases, the diameter of the rod remained unchanged, within the 
resolution of a high-quality micrometer. 

The second model system was an SCS-6 fiber (Textron Spe- 
cialty Materials, Lowell, MA) which has a —33 pan diameter 
pitch-based carbon core within a SiC sheath. Both the carbon 
core and the CVD SiC sheath have been well characterized 
using transmission electron microscopy. 20,21 The carbon core 
consists of a random mixture of small (1-50 nm) turbostatic 
carbon blocks. The SCS-6 fibers were exposed at two tempera- 
tures, 800° and 1000°C. At these temperatures, oxidation of the 
SiC sheath should be negligible. 22 Fiber lengths from 2.5 to 12 
cm, depending on the anticipated recession distance, were used. 
The furnace was heated rapidly to temperature and the speci- 
mens were held at the desired temperature for times to 50 h. 
For the longest fibers, the equilibrium temperature was — 10°C 
lower at the ends than the center. The samples were removed 
after cooling the furnace to room temperature. 

Following exposure, the recession lengths of the carbon cores 
were measured. The details of how this was accomplished are 
summarized elsewhere. 23 In short, the measurements were 
accomplished by mounting a section of the fiber on a plate with 
epoxy and grinding to approximately one-half of its original 
diameter. The recession length of the carbon core could then 
easily be measured using optical microscopy. Figure 5 shows 
typical optical micrographs of a ground fiber following 2 h of 
oxidation at 800°C. In both model systems, the receded carbon 
surface was planar following oxidation at high temperature or 
long times. However, at lower temperatures and shorter times, 
the carbon had a pointed morphology, as seen in the figure. In 
these cases, recession distances were recorded from the end of 
the sheath to the tip of the receded carbon. 




Fig. 5. Optical micrographs of sectioned and polished SCS-6 fiber 
oxidized for 2 h at 800°C. 


Solely for experimental convenience, the two systems were 
oxidized under different atmospheric conditions, the carbon 
rods in nominally static air and the SCS-6 in flowing (—6 cm/s) 
oxygen. In both, the atmosphere external to the pore is well 
mixed. In the case of nominally static air, this is due to thermal 
convection. The ambient pressure was the same in all experi- 
ments (1 atm). Therefore, the only difference of significance is 
the value of x, which is either 0.2 or 1 . 

Reproducibility of the data was checked in two ways. On the 
SCS-6 fibers, measurements from both ends of each oxidized 
fiber were taken and found to be consistent. In both model 
composite systems, selected experiments were repeated, and 
again a high degree of consistency was observed. 

IV. Results 

Predictions for recession distances in the two very different 
microcomposite systems were calculated using the linear- 
parabolic relationship of Eq. (1). It was assumed that the 
reaction was CO, and k p was calculated using Eq. (76) with 
appropriate values for d, x, T, and P. The linear rate constant is 
dependent on the physical and chemical characteristics of the 
carbon being oxidized; a different k x is required for each system, 
and can be obtained only by fitting experimental data. 

For the carbon rod, it was determined that k t was so large 
that the recession behavior was accurately modeled as purely 
parabolic. Figure 6 presents both the experimental data and 
curves calculated for the corresponding temperatures using 
Eq. (76). It can be seen that the model, which has no adjustable 
parameters, gives very accurate predictions. For this sys- 
tem, recession rate depends very weakly on temperature, as 
expected, since it is controlled by gas phase diffusion. 

As mentioned in Section III, the receding interface did not 
always remain planar, but developed a pointed morphology at 
lower temperatures. The presence or absence of the pointed tip 
does not appear to affect the agreement between the model and 
the experimental data. This is not surprising. The principal 
effect of such a tip shape is to expose additional surface area 
for reaction, effectively increasing k u which is already appar- 
ently so large that it can be neglected. 

Very different behavior was observed with the SCS-6 cores. 
Figure 7 shows a comparison of the experimental and predicted 
recession distances for the SCS-6 fibers at 800° and 1000°C. 
Again, the predicted recessions were calculated using Eq. (1) 
and obtaining k p from Eq. (76). In this case, however, assuming 
purely parabolic behavior resulted in a gross overestimation of 
recession distances. In other words, the oxidation rate was 
slower than would be expected based on diffusion control, and 
oxidation must be limited, over a significant amount of time, by 


Recession 

Distance 

(meters) 


1 hour 1 day 



Fig. 6. Comparison of predicted recession distances to experimental 
observations for 3.175 mm diameter carbon rods in alumina sheaths (A 
775°C, O 900°C, □ 1000°C, O 1200°C). 
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the chemical reaction at the receding carbon surface. The value 
of k x that gave the best fit to the 800°C data was 2.5 X 10 -6 m/s. 

V. Discussion 

Carbon oxidation is, of course, generally recognized to be 
fast, but the general analytical model allows quantitative predic- 
tion of recession rates and important variables to be analyzed. 
It is perhaps most illuminating to initially consider the transi- 
tion from linear to parabolic kinetics. Although it is recognized 
that the transition is broad (i.e., occurring over a wide range in 
time), it is possible to define the transition as the recession 
distance, or equivalently the oxidation time, at which the two 
terms on the left-hand side of Eq. (1) are equal. Equating these 
terms and substituting using Eqs. (6) and (59), i.e., assuming 0 2 
diffusion in and CO diffusion out, gives 

£,ra„s = * P /*i = (4/x) In (1 + X)(D/K r ) (77) 

Thus, the transition occurs at a depth which depends on the type 
of carbon (through K t ), the partial pressure of oxygen, the total 
system pressure, and temperature. The transition depth for air 
(X = 0.2) is ~26% of the depth for pure O z (x = 1). The 
transition occurs earlier (smaller values of (j trans ) at lower system 
pressure because k x is proportional to pressure, but k p (in the 
molecular regime) is independent of pressure. The temperature 
dependence is primarily through the exponential function, in 
K r , and the transition is delayed (occurs at larger values of £; lnms ) 
at lower temperatures. 

The fits to the experimental data allow transition distances 
to be estimated. For the case of the carbon rods, the transition 
distances are negligibly small at all temperatures investigated, 
on the order of 10 nm or less. During fitting it was observed 
that using k x s small enough to correspond to larger transition 
distances caused poorer fits. In contrast, the SCS-6 cores were 
best fit using values of k x that yield transition distances in excess 
of 5 cm at 800°C, and nearly 1 cm at 1000°C. Thus, it is entirely 
plausible for recession kinetics to fall into either of the two 
limiting cases of Eq. (1), for distances of engineering interest. 
It depends nearly entirely on the chemical and physical proper- 
ties of the carbon. 

The results of the analysis given in Section II can be summa- 
rized as follows. When the recession distance, is small, such 
that £ << £ trans , the recession rate will be linear* inversely 
proportional to pressure, exponentially dependent on tempera- 
ture, and the type of carbon will affect the oxidation rate. In the 
other limiting case, £ >> £ trans , recession will be parabolic, 


1 day 



Fig. 7. Comparison of predicted recession distances to experimental 
observations for 33 jam carbon core in SCS-6 SiC sheath. 


independent of the type of carbon, and weakly dependent on 
temperature and the nature of the gas phase reaction sequence. 
Lastly, the parabolic recession rate will be independent of pres- 
sure in the molecular regime, but directly proportional to pres- 
sure in the Knudsen regime. 

Knudsen effects are predicted to be very modest; i.e., carbon 
recession will be rapid even when the effective gas phase diffu- 
sion coefficient is D k . Consider first the behavior anticipated 
for carbon coatings with thicknesses in the range of 0. 1-1.0 
pm. Referring to Fig. 3, it can be seen that though both oxidize 
more slowly than that predicted using the molecular approxima- 
tion, the effect is small for both. The values of k p at 1000°C and 
1 atm for the 1.0 and 0.1 pm coatings are roughly 80% and 
30%, respectively, of the molecular prediction. Although these 
values are lower, they do not qualitatively change the situation. 
Substitution of the reduced values for k p into Eq. (77) still 
yields £ lrans on the order of centimeters. Furthermore, if the 
ambient pressure is raised to 15 atm, the 1.0 pm coating falls 
completely into the molecular regime and the 0. 1 pm coating is 
predicted to have a k p suppressed by less than 15%, relative 
to the molecular case. Thus, in heat engine applications of 
composites, Knudsen effects are essentially of no consequence. 

The fact that Knudsen diffusion is not slow has other 
important implications for composites. One class of alternate 
interface coatings is a multilayer approach in which concentric 
cylinders of carbon and silicon carbide are deposited. In such 
systems, the thickness of the carbon coatings is perhaps an 
order of magnitude thinner than simple coatings, on the order 
of tenths of micrometers rather than a few micrometers. One 
purported advantage of this approach is that the oxidation kinet- 
ics will be retarded. The above discussion makes it clear that 
reducing the thickness of the carbon to this range will not be 
sufficient to suppress oxidation to an acceptable rate. 

Also, it is generally recognized that oxidation under stress of 
CFCCs is much more rapid than static oxidation. There are two 
possible contributors to an acceleration of oxidation: increasing 
the matrix-crack opening displacement and oxidant ingress into 
the debond region at the fiber interface. 

Increasing the crack opening displacement is unlikely to be 
of much consequence for two independent reasons. First, the 
rate of recession will usually be controlled by the reaction 
at the receding carbon interface; thus, facilitating gas phase 
diffusion by opening the crack will not affect the rate. Second, 
even when recession is controlled by gas phase diffusion, the 
crack opening need only be modest, fractions of a micrometer, 
before the diffusion coefficient reaches the limiting value of the 
molecular regime. 

Oxygen ingress into the debond region represents an interest- 
ing situation. When a matrix crack passes through a composite, 
the interface between the intact fibers which span the crack and 
the matrix is debonded over a finite distance. If the debonded 
region has a finite opening, gas phase diffusion will rapidly lead 
to oxidation along the length of the debond (even if the opening 
leads to diffusion in the Knudsen regime). This leads to an 
instability, because oxidation will decouple the matrix from the 
fiber, preventing stress transfer, which indirectly results in a 
further length of interface being debonded. The newly 
debonded region is then subject to oxidation and this repeats in 
a continuous cycle. It is plausible that this mechanism causes 
the generally observed strong coupling between stress and 
oxidation. 


VI. Summary and Conclusions 

Oxidation of carbon in model composites is well described 
by a linear-parabolic law. The parabolic constant, k p , can be 
predicted using the kinetic theory of gases, whereas k x must be 
determined by fitting experimental data. The consequences of 
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this work for conditions of engineering interest (800°C < 
T < 1200°C and recession distances <1 mm) and experimental 
studies are many. Whether CO or C0 2 is formed directly, or 
whether C0 2 is the reaction product of sequential oxidation 
reactions, is an unimportant consideration for most applica- 
tions. For the case of a fiber-reinforced composite, it is unlikely 
that a carbon interphase can be fabricated thin enough to slow 
oxidation to a technologically acceptable rate. Once carbon is 
exposed to oxygen, via pores or matrix cracks, oxidation will 
proceed at a rapid rate, nearly independent of the magnitude of 
the pore/crack opening, the partial pressure of oxygen, or the 
total system pressure. The single, overriding, variable of conse- 
quence in moderating the kinetics is controlling the rate of 
chemical reaction at the receding carbon interface. 

APPENDIX 

List of Parameters Used in Model 


Symbol 

Definition 

Units 

C T 

Total concentration of gases 

mol/m 3 

c f 

Concentration of specie i at distance x 

mol/m 3 

c 0 M 

Concentration of oxygen in the 



environment 

mol/m 3 

c ox (4) 

Concentration of oxygen at carbon/ 



gas interface 

mol/m 3 

d 

Pore opening diameter 

m 

d' 

Annular pore opening thickness 

m 

D, 

Molecular diffusivity of specie i 

m 2 /s 

£> k 

Knudsen diffusivity for oxygen 

m 2 /s 

ji 

Flux of specie i 

mol/m 2 -s 

K r 

Reaction rate constant 

m/s 

K 

Linear recession rate constant 

m/s 

k Q 

Preexponent in Arrhenius form of k t 

m/s 

k p 

Parabolic rate constant for recession 

m 2 /s 


Molecular weight of specie i 

N/mol 

N 

Molar density of carbon 

mol/m 3 

P 

Pressure of the environment 

MPa 

Q 

Activation energy for carbon 



oxidation reaction 

J/mol 

r 

Rate of reaction 

mol/m 2 -s 

R 

Gas constant 

J/mol -K 

t 

Time 

s 

T 

Absolute temperature 

K 

V, 

Velocity of the specie i 

m/s 

e/k 

Energy of molecular interaction 



normalized by Boltzmann’s constant 

K 

°ij 

Collision diameter for species i,j 

m 

fl 

ij 

Collision integral (tabulated) for ij 


ji 

flux of specie i 

mol/m 2 -s 

X 

Distance from the surface into the pore 



(Fig. 1) 

m 

i 

Recession distance 

m 

X 

Fraction defining partial pressure of 
oxidant 
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